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TECHNICAL MEMORANDUM 


PATRAN-STAGS TRANSLATOR (PATSTAGS) 
INTRODUCTION 


PATSTAGS will translate PATRAN finite element model data into STAGS (Structural 
Analysis of General Shells) input records. PATSTAGS reads the data Irom a PATRAN neutral iile 
and writes STAGS input records into a STAGS input file and a UPRESS data file. The translator 
will not translate all PATRAN neutral file packets, nor will it write a complete STAGS input deck. 
The STAGS input deck must be edited after translation and the appropriate control records added 
as specified by the STAGS user manual. The following PATRAN neutral file packets are read and 
translated to the corresponding STAGS records. 


PATRAN neutral file packet 

To 

STAGS input record 

25 - Title Card 



26 - Summary Data 



1 - Node Data 


S-l User Point Record 

S-2 Auxiliary Coordinate Record 

2 - Element Data 


T-2 Beam Element Record 
T-3 Triangular Element Record 
T-4 Quadrilateral Element Record 

6 - Distributed Loads 


UPRESS input record 

7 - Node Forces 


U-3 Load Definition Record 

8 - Node Displacements 


S-l User Point Record 


FILES 


PATSTAGS uses three files: the PATRAN neutral file to be translated, a STAGS input file, 
and a STAGS pressure data file. PATSTAGS will prompt the user for the name of the neutral file 
to be translated and the desired names of the STAGS files to be created. The file names may be 
up to 40 characters in length. The STAGS input file created will contain the STAGS S-l, S-2, 

T-2, T-3, T-4, and U-3 input records. The pressure data file created will contain the element live 
pressure data used by the STAGS subroutine UPRESS. 




NODAL DATA 


Nodal data is read from the PATRAN neutral file packets 1 (Node Data) and 8 (Node Dis- 
placements). The data is written as S-l cards in the STAGS input file. S-2 cards are also created if 
cylindrical coordinate system 1 is used for the node coordinates and the node constraints in the 
PATRAN model. Up to 5,000 nodes may be translated. To translate a larger model, the array 
dimensions will need to be increased in PATSTAGS.FOR and the program relinked. 

The translator has the capability to define an auxiliary cylindrical coordinate system on a 
S-2 card for use in applying boundary constraints. The following requirements must be met before 
an S-2 card is written. 

1. The node and nodal constraints must both be defined in PATRAN in cylindrical 
coordinate system 1 . 

2. The global X axis must coincide with the cylindrical Z axis, i.e., the model must be a 
shell of revolution about the global X axis. 

The auxiliary coordinate system for each node is defined on S-2 cards with two points. The 
points are defined as follows: 


1 . The node radius is calculated as 


Rad n — V K - +Z„ 2 


where Y„ and Z„ are the global y and z coordinates of the node. 
2. Point 1 coordinates are calculated as 


X, = 0.0 


Yi 


Y„ 

Rad n 


Z, = 


Rad„ 


where subscript 1 refers to the cylindrical coordinate system point 1 , and subscript n refers to the 
node. 
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3. Point 2 coordinates are calculated as 


X -> = 0.0 


Y-> = - 


Rad,, 


Z , = 


Rad, , 


where subscript 2 refers to cylindrical coordinate system point 2 and subscript n refers to the node. 

Figure 1 shows the coordinate systems attached to four nodes in the different quadrants. 

The coding for the auxiliary coordinate system is contained in the 9200 block of the 
PATSTAGS FORTRAN file. This can be easily changed if a different axis of revolution is desired 
for the model. 

Nodal constraints only are supported. Specified nodal displacements or rotations which 
would require a U-3 record to be written are not currently supported. 



Figure I . Auxiliary coordinate systems. 
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ELEMENT DATA 


Element data is read from the PATRAN neutral file packet 2 (Element Data). The data is 
written as T-2, T-3, and T-4 cards in the STAGS input file. Bar, triangular, and quadrilateral 
elements are supported. Up to 5,000 elements may be translated. To translate a larger model, the 
array dimensions in PATSTAGS.FOR must be increased and the program relinked. 

The program will interactively write the number of elements which have been read and will 
prompt the user for the desired values of ILIN, which governs geometic nonlinearity, and IPLAS, 
which governs material nonlinearity. The entered response will be used in all T-2, T-3, and T-4 
records written. The property ID numbers entered in PATRAN will be used for 1CROSS in the T-2 
card and IWALL in the T-3 and T-4 cards. XSI, EC4, and ECQ in the T-2 records, ZETA and 
ECZ in the T-3 records, and ZETA, ECZ, INTEG, and 1PENL in the T-4 records are all set to 
zero. Any of these defaults can be changed by editing the appropriate format card in block 9400 of 
the PATSTAGS.FOR file and relinking the program. The program will next interactively write the 
number of beam elements read and prompt the user for the desired beam element code number. 

The triangular and quadrilateral elements are handled similarly. 

At the end of each element record in the STAGS input file, the record type, PATRAN 
element number, and STAGS element number are written for the user’s information. 

The X-Y plane of all bar elements must be defined using the node option in PATRAN. The 
vector and grid option are not supported. 


FORCE DATA 


Force data is read from the PATRAN neutral file packet 7 (Node Forces). The data is 
written as U-3 cards in the STAGS input file. 

The translator will support only one force component per node, per load set. If more than 
one force component is needed on a node, they should be defined in PATRAN as belonging to dif- 
ferent load sets (i.e., the X component in load set I , the Y component in load set 2, the Z compo- 
nent in load set 3, etc.). 


PRESSURE DATA 


Pressure data is read from the PATRAN neutral file packet 6 (Distributed Loads). The data 
is written into a pressure data file, which is then used with the UPRESS subroutine. The translator 
interactively prompts the user for the desired name of the pressure data file. 
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The translator will support only one pressure component per element, per load set. If more 
than one pressure component is needed on an element, they should be defined in PATRAN as 
belonging to different load sets (i.e. , the X component in load set 1, the Y component in load set 
2, the Z component in load set 3, etc.). 

This data file is formatted to be used with the subroutine UPRESS listed in appendix A. 
This subroutine will need to be linked to STAGS before running the analysis. The pressure data 
file is called by the subroutine UPRESS as unit 17, so an assign statement is needed to assign the 
pressure data file to unit 17 before running the analysis. 
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APPENDIX A 


SUBROUTINE UPRESS 
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SUBROUTINE UPRESS(T, PA, PB,IUNIT,IELT,X,Y.Z,LIVE,PRESS) 


C 

C THIS SUBROUTINE, WHEN LINKED WITH STAGS, WILL READ LIVE 
C PRESSURE DATA FROM THE FORMATED PRESSURE FILE CREATED 
C BY PATSTAGS. 

C 

C WRITTEN BY NEIL OTTE 

C MARSHALL SPACE FLIGHT CENTER 

C ED-24 STRUCTURAL STRENGTH BRANCH 

C (205) 544-7231 

C 

DIMENSION STID(5000),EPRESS(5000) 

INTEGER STID 
IF (M.EQ.O)THEN 
READ (17,10)NPRESS 

10 FORMAT (15) 

DO 100 I=1,NPRESS 

READ (17,11) STID(I),EPRESS(I) 

11 FORMAT (I5.F10.4) 

100 CONTINUE 

WRITE (6,12) 

12 FORMAT (IX, ’SUBROUTINE UPRESS WRITEN BY NEIL OTTE',/, IX, 

* HAS BEEN USED’) 

ELSE 

CONTINUE 
END IF 
K = 1 

200 CONTINUE 

IF (STID(K).EQ.IELT)THEN 
PRESS = EPRESS 
LIVE = 1 
ELSE 
K = K+l 
GO TO 200 
ENDIF 
M = 99 
RETURN 
END 
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APPENDIX B 
PATSTAGS.FOR 
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OTocxjoooooooccocmm^ 

c 

c 

C PATRAN - STAGS TRANSLATOR 

C WRITTEN BY: NEIL OTTE 

C STRUCTURAL STRENGTH BRANCH ED-24 

C MARSHALL SPACE FLIGHT CENTER 

C (205) 544-7231 

C 

C THIS PROGRAM WILL READ A PATRAN NEUTRAL FILE AND CREATE A STAGS INPUT 

C DECK, AS WELL AS A LIVE PRESSURE DATA FILE TO BE READ BY THE UPRESS 

C SUBROUTINE. THIS PROGRAM WILL NOT TRANSLATE ALL PATRAN NEUTRAL 

C FILE PACKETS NOR WILL IT WRITE A COMPLETE STAGS INPUT DECK. SEE THE 

C USERS MANUAL FOR MORE INFORMATION. 

C 

c 

oooocriiximxxiixxxxDcccxxixiixmDcmxmxiixmTr block 9000 ccccaxcooocc 

c 

C DIMENSION AND INITIALIZATION BLOCK 

C 

DIMENSION NID(5000),X(5000),Y(5000),Z(5000),CID(5000),XD(5000), 

* YD(5000),ZD(5000),UX(5000),UY(5000),UZ(5000),NDCID(5000), 

* NDID(5000),EBPID(5000),EBLNODE1(5000),EBLNODE2(5000), 

* EBNODER(5000),ETPID(5000),ETLNODE1(5000),ETLNODE2(5000), 

* ETLNODE3(5000),EQPID(5000)£QLNODE1(5(X)0),EQLNODE2(5000), 

* EQLNODE3(5000),EQLNODE4(5000),BID(5000),TID(5000), 

* QID(5000),STID(5000),LID(5000),PDATA(5000),LE(5000), 

* LD(5000),FN(5000),FCID(5000),FD(5000),FDATA(5000), 

* FLAX(5000) 

INTEGER PT,EID,EIV,EKC,EN1 ,EN2,EBNODES,EBCONFIG,EBPID,EBCEID, 

* EBLNODE1 ,EBLNODE2£BNODERTTNODESJETCONFIGTTPID,ETCEID, 

* ETLNODE1 JTIUNODE2£TLNODE33QNODESJEQCONHG£QPID£QCEID, 

* EQLNODE1 EQLNODE2,EQLNODE3£QLNODE4,BID,TID,QID,STID,CID, 

* CONFIG, PSPC1,PSPC2,PSPC3,PSPC4,PSPC5,PSPC6,XD, YD ,ZD, 

* UX,UY,UZ,FD,FN,FCID,FICOMPl,FICOMP2,FICOMP3,FICOMP4, 

* FICOMP5,FICOMP6,FE,FLAX 

CHARACTER GTYPE*10,INFILE*40,OUTFILE*40,TITLE*80,DATE*12,TIME*8, 

* VERS ION* 12, PRESS URE*40 
C 

C INITIALIZATION 

C 

IB = 0 
IT = 0 
IQ = 0 
C 

caDocooaxraxxraxraxxxcctt^ block 9 ioo coooccccccooc 

c 

C FILE SETUP BLOCK 

C 

C PROMPT FOR INPUT AND OUTPUT FILE NAMES 

C 

WRITE (5,9100) 

9100 FORMAT (2X, 'ENTER THE NAME OF THE NEUTRAL FILE TO BE TRANSLATED ) 

READ (6,9110) INFILE 
9110 FORMAT (A40) 


PRECEDING PAGE BLANK NOT 

PRECEDING page BUNK not 


13 

INI£NTiaiiAU.K BLANK 



WRITE (5,9101) 

9101 FORMAT (2X, 'ENTER THE DESIRED NAME OF THE STAGS INPUT FILE') 

READ (6,9110) OUTFILE 

WRITE (5,9102) 

9102 FORMAT (2X, ENTER THE DESIRED NAME OF THE STAGS PRESSURE FILE ) 

READ (6,9110) PRESSURE 

C 

C OPEN THE INPUT AND OUTPUT FILES 

C 

OPEN (UNH^.FILE^NFILE.STATUS^ UNKNOWN') 

OPEN (UNIT =8 ,FILE=OUTFILE ,ST ATUS=’UNKN OWN') 

OPEN (UNIT= 1 0,FTLE=PRESSURE, STATUS-UNKNOWN') 

C 

aXXXXXXXXXXXOeOCXXXXXXXXXYXXXXYXXXXXXXXXXTXXXXXTXi: BLOCK 250 oococxxxxxxxxr 

c 

c TITLE CARD BLOCK 

c 

C READ PACKET TYPE 25 - TITLE CARD 

C 

READ (7,250) PT,ID,IV,KC, TITLE 
250 FORMAT (I2,3I8,/,A80) 

C 

cxxxxxxxxxxxraxxxxxxxxxixxxxxxrxxxixxxxxTxxxxxxxxxx: BLOCK 260 OOOCOOCCCCCCC 
c 

C SUMMARY DATA BLOCK 

c 

C READ PACKET TYPE 26 - SUMMARY DATA 

C 

READ (7,260) PT.ID, I V,KC,N1,N2,N3,N4,N5,DATE,TIME, VERSION 
260 FORMAT (12,818,/, 12A.8A.12A) 

C 

(XXXXXXXXXXXXYXXXXTXXXXXXXLXXIXXXXXXXXXTXXXXXXXXXXXC block io axxxxccxxxccc 
c 

c NODE COORDINATE BLOCK 

C 

C READ PACKET TYPE 0 1 - NODE COORD. DATA 

C 

DO 100 1=1, N1 

READ (7,10) PT,NID(I),NIV,NKC,X(I),Y(I),Z(I),ICF,GTYPE,NDF, 

* CONFIG, CID(I),PSPC1,PSPC2, PS PC3.PSPC4, PSPC5.PS PC6 

10 FORMAT (I2,3I8,/,3E16.9,/,I1 ,1 A1 ,3I8,2X,6I1) 

100 CONTINUE 
C 

aXXXXtXXaXXXIXXXXXXXmXXXXXXXXXXXXXXXXXXXXXXXXIXE BLOCKS (XOOOCOOOOOOOC 

c 

C ELEMENT DATA BLOCK 

c 

C READ PACKET TYPE 02 - ELEMENT DATA 

C 

DO 200 1=1, N2 

READ (7,20) PT,EID,EI V,EKC,EN 1 ,EN2 
20 FORMAT (12,818) 

C 

C READ BEAM ELEMENT DATA 

C 
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IF (EIV.EQ.2)THEN 
IB = IB+1 
BID(IB) = EID 

READ (7,21) EBN0DES,EBC0NFIG,EBPID(1B),EBCEID,EBTH1.EBTH2, 

* EBTH3,EBLNODEl(IB),EBLNODE2(IB) 

21 FORMAT (4I8.3E16.9A2I8) 

EBNODER(IB) = EN2 

ELSE 

CONTINUE 
END IF 
C 

C READ TRIANGULAR ELEMENT DATA 

C 

IF (EIV.EQ.3)THEN 
IT = IT+1 
TID(IT) = EID 

READ (7,22) ETNO DES ,ET CONFIG ,ETPID(IT),ETCEID ,ETTH 1 JETT H2 .ETTH3 , 

* ETLNODEl(IT),ETLNODE2(IT),ETLNODE3(IT) 

22 FORMAT (4I8.3E 1 6.9./.3I8) 

ELSE 

CONTINUE 
END IF 
C 

C READ QUADRALATERIAL ELEMENT DATA 

C 

IF (EIV.EQ.4)THEN 
IQ = IQ+I 

QID(IQ) = EID 

READ (7,23) EQNODES,EQCONFIG,EQPID(IQ) JEQCEIDJEQTH 1 .EQTH2.EQTH3, 

* EQLNODE1(IQ),EQLNODE2(IQ),EQLNODE30Q),EQLNODE4{IQ) 

23 FORMAT (4I8,3E16.9,/,4I8) 

ELSE 

CONTINUE 

ENDIF 

200 CONTINUE 
C 

C CREATE MATRIX OF STAGS ELEMENT ID VS. PATRAN ELEMENT ID 

C 

IF (IB.GT.O)TIIEN 
DO 201 L =1,IB 
STID(L) = BID(L) 

201 CONTINUE 
ELSE 

CONTINUE 

ENDIF 

IF (IT.GT.O)THEN 
M = IB+1 
N = IB+IT 
DO 202 L =M,N 
K = L-IB 

STID(L) = TID(K) 

202 CONTINUE 
ELSE 

CONTINUE 

ENDIF 
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IF (IQ.GT.O)THEN 
M = IT+IB+1 
DO 203 L =M,N2 
K = L-IB-IT 
STID(L) = QID(K) 

203 CONTINUE 
ELSE 

CONTINUE 

ENDIF 

C 

ctxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx: block 30 ccaxxxxxxxxxxxxxr 
c 

C SKIP BLOCK 

c 

C SKIP PACKET TYPES 3, 4, AND 5 

c 

300 CONTINUE 

READ (7,30) PT,ID,IV,KC 
30 FORMAT (12,818) 

IF (PT.EQ.3) THEN 
DO 301 1=1, KC 
READ (7,*) 

301 CONTINUE 
GOTO 300 

ELSE 

CONTINUE 

ENDIF 

IF (PT.EQ.4) THEN 
DO 400 I =1,KC 
READ (7,*) 

400 CONTINUE 
GOTO 300 
ELSE 

CONTINUE 

ENDIF 

IF (PT.EQ.5) THEN 
DO 500 1=1, KC 
READ (7,*) 

500 CONTINUE 
GOTO 300 
ELSE 

CONTINUE 

ENDIF 

IF (PT.EQ.6)THEN 
C 

OOCCCXXXXXXXXXXXXXXXXXXXXXXXXXJXXXXXXXJJXXXXXXXX: BLOCK 60 OOOCXXDOCmXXDOOCC 

c 

C PRESSURE BLOCK 

c 

C READ PACKET TYPE 6 - PRESSURE LOAD 

C 

K = 1 

600 CONTINUE 
LID(K) = ID 

READ (7,61) LTYPE,LEFLAG,LGFLAG,LICOMPl,LICOMP2,LICOMP3,LICOMP4, 


* LICOMP5,LICOMP6,LNODEl JLNODE2,LNODE3,LNODE4,LNODE5,LNODE6, 

* LNODE7,LNODE8,LNFE 

61 FORMAT (1711,12) 

READ (7,62) PDATA(K) 

62 FORMAT (5E16.9) 

C 

C FIND STAGS ELEMENT ID FOR THIS LOAD 

C 

1 = 0 

601 CONTINUE 

I = 1+1 

IF (STID(I).EQ.LID(K))THEN 
LE(K) = I 
ELSE 
GO TO 601 
ENDIF 
C 

C FIND PROPER LOAD DIRECTION 

C 

IF (LICOMP1 .EQ. 1 )LD(K)= 1 
IF (LICOMP2 .EQ . 1 )LD(K)=2 
IF (LICOMP3.EQ.l)LD(K)=3 
IF (LICOMP4.EQ. 1 )LD(K)=4 
IF (LICOMP5.EQ.l)LD(K)=5 
IF (LICOMP6.EQ. 1 )LD(K)=6 
KMAX = K 
C 

C READ NEXT HEADER 

C 

READ (7,60) PT,ID,IV,KC 
60 FORMAT (12,818) 

IF (PT.EQ.6) THEN 
K = K+l 
GO TO 600 
ELSE 

CONTINUE 

ENDIF 

ELSE 

CONTINUE 

ENDIF 

C 

occxDcccmDooocaDOoaxxxTxxxxxxarixxrxixxixmxc blockx ooccccccccoocaxxx: 
c 

c FORCE BLOCK 

C 

C READ PACKET TYPE 7 - FORCES 

C 

IF (PT.EQ.7)THEN 
M = 1 

700 CONTINUE 
FN(M) = ID 

READ (7,70) FCID(M),FICOMPl,FICOMP2,FICOMP3,FICOMP4,FICOMP5, 

* FICOMP6 
70 FORMAT (18,611) 

READ (7,71) FDATA(M) 


71 FORMAT (El 6.9) 

C 

C SET UP THE PROPER COORDINATE NUMBER 

C 

IF (FCID(M).EQ.l) THEN 
FLAX(M) = 0 
ELSE 

FLAX(M) = 1 
END IF 
C 

C FIND THE PROPER LOAD DIRECTION 

C 

IF (FICOMP1.EQ.1) FD(M)=1 
IF (FICOMP2.EQ.1) FD(M)=2 
IF (FICOMP3.EQ.1) FD(M)=3 
IF (FICOMP4.EQ.1) FD(M)=4 
IF (F1COMP5.EQ.1) FD(M)=5 
IF (F1COMP6.EQ.1) FD(M)=6 
MMAX = M 
C 

C READ NEXT HEADER 

C 

READ (7,72) PT,ID,IV,KC 

72 FORMAT (12,818) 

IF (PT.EQ.7) THEN 

M = M+l 
GO TO 700 
ELSE 

CONTINUE 

ENDIF 

ELSE 

CONTINUE 

ENDIF 

IF (PT.EQ.8)THEN 
C 

oooocoxocoocixiTxiiixxxxmDcixixxximiixmiix: blockso cooccDoooaDoooococ 
c 

C NODE CONSTRAINT BLOCK 

C 

C READ PACKET TYPE 8 - NODE CONSTRAINTS 

C 

800 CONTINUE 
I = ID 

NDID(I) = I 
IF (KC.EQ.3)THEN 

READ (7,81) NDCID(I),XD(I),YD(I),ZD(I),UX(I),UY(I),UZ(I), 

* DDATA1 ,DDATA2,DDATA3,DDATA4,DDATA5,DDATA6 

81 FORMAT (18,61 1 ,/,5E 1 6.9./.5E 1 6.9) 

ELSE 

READ (7,82) NDCID(I),XD(I),YD(I),ZD(I),UX(I),UY(I),UZ(I), 

* DDATA1 ,DDATA2,DDATA3,DDATA4,DDATA5 

82 FORMAT (I8,6I1,/,5E16.9) 

END IF 

READ (7,80) PT,ID,IV,KC 
80 FORMAT (12,818) 
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IF (PT.EQ.8) THEN 
GO TO 800 
ELSE 

CONTINUE 

ENDIF 

ELSE 

CONTINUE 

ENDIF 

C 

cxxxLCLTxuccLmxccmxTXXirxxxcxxxmxtxxrcoQC block 9200 cocoocoooccoooooc 

NODE DATA WRITE BLOCK 

WRITE NODAL DATA IN STAGS FORMAT 

DO 9200 I=1,N1 
IF (XD(I).EQ.0)THEN 
XD(I) = 1 
ELSE 
XD(I) = 0 
ENDIF 

IF (YD(I).EQ.0)THEN 
YD(I) = 1 
ELSE 
YD(I) = 0 
ENDIF 

IF (ZD(I).EQ.0)THEN 
ZD(I) = 1 
ELSE 
ZD(I) = 0 
ENDIF 

IF (UX(I).EQ.0)THEN 
UX(I) = 1 
ELSE 
UX(I) = 0 
ENDIF 

IF (U Y (I).EQ.0)THEN 
UY(I) = 1 
ELSE 
UY(I) = 0 
ENDIF 

IF (UZ(I).EQ.0)THEN 
UZ(I) = 1 
ELSE 
UZ(I) = 0 
ENDIF 

WRITE (8,9201) NIDCO.O.O.O.XCU.YCU.ZCO.XDCIJ.YDCO.ZDCI), 
UX(I),UY(I),UZ(I),CID(I) 

9201 FORMAT (IS.lXJLlX.Il.lX.IUX^AlXJ^AlX.FgAlX, 

* 3I1,1X,3I1,1X,I1,10X;$ S-l’) 

IF (CID(I).EQ. 1)THEN 
RAD = SQRT(Y(I)**2+Z(I)**2) 

IF (RAD.EQ.O)THEN 
XAY = 1.0 
XAZ = 0.0 


19 



ELSE 

XAY = (Y(I)/RAD) 

XAZ = (Z(I)/RAD) 

YAY = -(Z(I)/RAD) 

YAZ = (Y (I)/RAD) 

END IF 

WRITE (8,9202) 0.0,XAY,XAZ,0.0,YAY,YAZ 
9202 FORMAT (1X,F9.4,',',F9.4,7F9.4,1X,F9.4,\',F9.4, 

* ',',F9.4,5X,’$ S-2') 

ELSE 

CONTINUE 
END IF 

9200 CONTINUE 
C 

cooccccccojoocooooooooccooooooooooocooooooooc block 9300 cxxxxxxxxxmxxxooc 

c 

C INTERACTIVE PROMPTS BLOCK 

C 

C INTERACTIVE PROMPTS 

C 

C PROMPT FOR ANALYSIS TYPE 

C 

WRITE (6,9300) N2 

9300 FORMAT (2X. DATA FOR',I7,2X, 'ELEMENTS HAVE BEEN READ.7,/,2X, 

* PLEASE ENTER 0 FOR NONLINEAR STRAIN-DISPLACEMENT RELATIONS', 

* /,2X,'OR 1 FOR LINEAR STRAIN-DISPLACEMENT RELATIONS.') 

READ (5,9310) ILIN 

9310 FORMAT (110) 

WRITE (6,9301) 

9301 FORMAT (2X, PLEASE ENTER 0 FOR ELASTIC BEHAVIOR OR 1 FOR7,2X, 

* ’PLASTICITY EFFECTS') 

READ (5,9310)IPLAS 

C 

C PROMPT FOR BEAM ELEMENT TYPE 

C 

IF (IB.GT.O)THEN 
WRITE (6,9302) IB 

9302 FORMAT (2X/DATA FOR’,I7,2X,'BEAM ELEMENTS HAVE BEEN READ.7.2X, 

* PLEASE ENTER THE DESIRED BEAM ELEMENT CODE NUMBER.') 

READ (5,9310) KBM 

ELSE 

CONTINUE 

ENDIF 

C 

C PROMPT FOR TRIANGULAR BEAM ELEMENT TYPE 

C 

IF (IT.GT.O)THEN 
WRITE (6,9303) IT 

9303 FORMAT (2X,'DATA FOR’,I7,2X, 

* TRIANGULAR ELEMENTS HAVE BEEN READ.7.2X, 

* 'PLEASE ENTER THE DESIRED TRIANGULAR ELEMENT CODE NUMBER.') 

READ (5,9310) KTRI 

ELSE 


CONTINUE 

ENDIF 

C 

C PROMPT FOR QUAD BEAM ELEMENT TYPE 

C 

IF (IQ.GT.O)THEN 
WRITE (6,9304) IQ 

9304 FORMAT (2X,'DATA FOR', 17, 2X, QUAD ELEMENTS HAVE BEEN READ.7,2X, 

* 'PLEASE ENTER THE DESIRED QUAD ELEMENT CODE NUMBER. ) 

READ (5,9310)KQUAD 

ELSE 

CONTINUE 

ENDIF 

C 

coocccccoxciiixmDCtixxxixixiixixmxmxixrc block smoo ooocxxdcocccccccooc 

c 

c ELEMENT DATA WRITE BLOCK 

C 

C WRITE ELEMENT DATA IN STAGS FORMAT 

C 

C WRITE BEAM ELEMENT DATA 

C 

IF (IB.GT.O)THEN 
DO 9400 1=1, IB 

WRITE (8,9401) EBLNODEl(I),EBLNODE2(I),EBNODER(I),KBM,EBPID(I), 

* ILIN,IPLAS,BID(I),I 

9401 FORMAT ( 1 X,I4, 1 X,I4, 1 X,I4, 1 X,I4, 1 X,I3, 1X,'0.\1 X.'O.'.lX.’O.'.l X, 

* I1,1X,I1,10X,'$ T-2',1X,I5,3X,I5) 

9400 CONTINUE 

ELSE 

CONTINUE 

ENDIF 

C 

C WRITE TRIANGULAR ELEMENT DATA 

C 

IF (IT.GT.O)THEN 
DO 9402 1=1, IT 
M = I+IB 

WRITE (8,9403) ETLNODEl(I),ETLNODE2(I),ETLNODE3(I),KTRI,ETPID(I), 

* ILIN,IPLAS,T1D(I),M 

9403 FORMAT (1X,I4,1X,I4,1X,I4,1X,I4,1X,I3,1X,'0.',1X,’0.’, 

* 1X,I1,1X,I1,10X,'$ T-3',1X,I5,3X,I5) 

9402 CONTINUE 
ELSE 

CONTINUE 

ENDIF 

C 

C WRITE QUAD ELEMENT DATA 

C 

IF (IQ.GT.O)THEN 
DO 9404 1=1, IQ 
M = I+IB+IT 

WRITE (8,9405) EQLNODEl(I),EQLNODE2(I),EQLNODE3(I),EQLNODE4(I), 

* KQUAD,EQPID(I),1LIN,IPLAS,QID(I),M 

9405 FORMAT (lX.RlX.RlXRlX.RlX.RlX.B.lX.’O.'.lX.'O.'.lX, 
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* I1,1X,I1,1X,'0',1X,'0M0X,'$ T-4’,1X,I5,3X,I5) 

9404 CONTINUE 

ELSE 

CONTINUE 
END IF 
C 

cocxraxxxxxraDcxxmxc^^ block9soo (XDXxxxmxmxrrc 

c 

C FORCE WRITE BLOCK 

C 

C WRITE FORCES TO STAGS FORMAT 

C 

DO 9500 M=1,MMAX 

WRITE (8,9501) FDATA(M),FD(M),FN(M),FLAX(M) 

9501 FORMAT (F10.3,1X,T,1X,I1,1X,I5,1X, , 0MX > 

* I1,5X,'$ U-3 FORCES’) 

9500 CONTINUE 

C 

ODCXXOOOCCOCXrXDCXDCXXXXXXXIXXXXXXXXXXXXXOCOOaDCXXI^ BLOCK 9600 COOCCCCCCOCC 

c 

C PRESSURE WRITE BLOCK 

C 

C WRITE PRESSURE DATA TO UPRESS FORMAT 

C 

WRITE (10,9600) KMAX 

9600 FORMAT (15) 

DO 9601 K=1,KMAX 

WRITE (10,9602) LE(K),PDATA(K) 

9602 FORMAT (15, FI 0.4) 

9601 CONTINUE 
STOP 
END 

C 

cooooooooooaooocxrraDoooooooocrxxccraDGOcxxrnDoocxmiooooooooociDocxDooooooooooGoc 
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